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Abstract. Sand pile models are dynamical systems emphasizing the phenomenon of Self 
Organized Criticality (SOC). From A'^ stacked grains, iterating evolution rules leads to some 
critical configuration where a small disturbance has deep consequences on the system, in- 
volving numerous steps of grain fall. Physicists L. Kadanoff et al inspire KSPM, a model 
presenting a sharp SOC behavior, extending the well known Sand Pile Model. In KSPM with 
parameter D we start from a pile of A*' stacked grains and apply the rule: D — 1 grains can 
fall from column i onto the D — 1 adjacent columns to the right if the difference of height 
between columns i and i+1 is greater or equal to D. We propose an iterative study of KSPM 
evolution where one single grain addition is repeated on a heap of sand. The sequence of grain 
falls following a single grain addition is called an avalanche. From a certain column precisely 
studied for D = 3, we provide a plain process describing avalanches. 

Keywords: Discrete dynamical system, self-organized criticality, sand pile model. 



1 Introduction 

Sand pile models were introduced in as systems presenting a critical self-organized 
behavior, a property of dynamical systems having critical points as attractors. In the scope 
of sand piles, starting from an initial configuration of stacked grains the local evolution 
of particles is described by one or more iteration rules. Successive applications of such rules 
alter the configuration until it reaches an attractor, namely a stable state from which no 
rule can be applied. SOC property means those attractors are critical in the sense that a 
small perturbation — adding some more grains — involves an arbitrary deep reorganization 
of the system. Sand pile models were well studied in recent years ([5]. [5]. [6] .115)). 



1.1 KadanofF sand pile model 

In |12] . Kadanoff proposed a generalization of classical models closer to physical behavior 
of sand piles in which more than one grain can fall from a column during one iteration. 
Informally, Kadanoff sand pile model with parameter D and N grains is a discrete dynam- 
ical system, which initial configuration is composed of stacks grains, moving in discrete 
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space and time according to a transition rule : if the height difference between column i 
and i + 1 is greater or equal to D, then D — 1 grains can fall from column i to the D — 1 
adjacent columns on the right (see figured]). 



> D 




Fig. 1: KSPM(D) transition rule. 



Sand pile models are specializations of Chip Firing Games (CFG). A CFG isplayed on 
a graph in which each vertex v has a load l{v) and a threshold t{v) = deg~^{vj^, and the 
transition rule is: if l{v) > t{v) then v gives one unit to each of its neighbors (we say v is 
fired). As a consequence, we inherit all properties of CFGs. 

Kadanoff sand pile is referred to a linear chip firing game in [11]. The authors show 
that the set of reachable configurations endowed with the order induced by the successor 
relation has a lattice structure, in particular it has a unique fixed point. Since the model 
is non-deterministic, they also prove strong convergence i.e. the number of iterations to 
reach the fixed point is the same whatever the evolution strategy is. The morphism from 
KSPM(3) to CFG is depicted on figure EJ 




Fig. 2: The initial configuration a of KSPM(3) is presented as a CFG where each vertex 
corresponds to a column (except the sink) seen as a difference of height. 



More formally, sand pile models we consider are defined on the space of ultimately null 
decreasing integer sequences. Each integer represents a column of stacked sand grains and 
transition rules describe how grains can move from columns. Let h = {Hq, hi, /12, . . . ) denote 
a configuration of the model, where each integer hi is the number of grains on column i. 
Configurations can also be given as height difference a = (cto, fxi, (T2, . . . ), where for all 
^ > 0, ai = hi — /ij+i. We will use this latter representation throughout the paper, within 
the space of ultimately null non-negative integer sequences. 

^ deg'^{v) denotes the out-degree of v. 
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Definition 1. The Kadanoff sand pile model with parameter D, KSPM(D), is defined by: 

— A set of configurations, consisting in ultimately null non-negative integer sequences. 

— A set of transition rules ; we have a transition from a configuration a to a configuration 
a' on column i, and we note a ^ a' if 

• (7[_-^ = (7,^1 + D-1 (fori^O) 

• <y'i = (^i — D , 

• '^'i+D-1 = + 1 

• c^'j = CTj for j ^ {i-l,i,i + D -1}. 

Remark that according to the definition of the transition rules, a condition for a' to be 
a configuration is that fij > D. 



1.2 Strategies and avalanches 

A basic property of the KSPM model is the diamond property. If there exists two distinct 
integers i and j such that a ^ a' and a ^ a'\ then there exists a configuration a'" such 
that a' a'" and a" A a'" . We note o" — )• a' when there exists an integer i such that 
a -V a'. We define the transitive closure A of — >•, and say that a' is reachable from a when 
a a . 

A strategy is a sequence s = (si, . . . , st). We say that a' is reached from a via s when 
(7 — ?• cr" A . . . fj' and we note a a' . We also say, for each integer t such that 
< t < r, that the column sj is fired at time t in s. (informally, the index of the sequence 
is interpreted as time). 

For any strategy s and any nonnegative integer i, we state \s\i = #{t|sj = i}. Let s^, 

s^ be two strategies such that a ^ and a A . We have the equivalence: [V i, \s^\i = 
\s^\i\ fj'^ = 0"^. A strategy s such that o" A cr' is called leftmost if it is the minimal 
strategy from a to a' according to lexicographic order. A leftmost strategy is such that at 
each iteration, the leftmost possible transition is performed. 

We say that a configuration a is stable, or a fixed point if no transition is possible 
from (T. As a consequence of the diamond property, one can easily check that, for each 
configuration cr, there exists a unique stable configuration, denoted by 7r((T), such that 
a A 7r(cT). Moreover, for any configuration a' such that a A a', we have it{(j') = vr(fT) (see 
[m for details). 

In this paper, we are interested in the iterative process defined below. Starting with 
no grain, we successively add a single grain on column 0, and make all the possible firings 
until a fixed point is reached. We denote by 7r(A;) the configuration obtained with this 
process using k grains (from the structure of KSPM described above, one easily checks 
that 7r(A:) = 7r((A:,0^)). 

Let (T be a configuration, a^^ is the configuration obtained by adding one grain on 
column 0. In other words, if o" = (ctq, ui, . . . ), then a^^ = (cjo + 1, (Ti, . . . ). 
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Formally the process is defined by 7r(0) = 0"^ and the recurrence formula: 

7r{TT{k - =TT{k). 

The k^^ avalanche s*^ is the leftmost strategy from 7r(/c — l)-'-" to n^k). The goal of 
the present paper is the description of avalanches. Informally, we want to describe what 
happens when a new grain is added in a previously stabilized sand pile. 

For D = 2, i.e. the classical SPM, this description is easy: the added grain moves 
rightwards until it arrives in a plateau. But, for D > 2, the situation is not so simple. We 
now state our results. 

— In the general case, we prove (Section [2]) the following properties: 

• Each column is fired at most once, 

• For any avalanche, as soon as an interval {L,L + 1,....,L + D — 1} of successive fired 
columns exists, the execution of the avalanche on the right part of this interval can 
be turned into a pseudo local and elementary process. 

Informally, that means that the knowledge of such an interval guarantees a regular 
behavior of the avalanche on its right part. 

— In the case when D = 3, we prove (Section [S]) the property below: 

• For each avalanche s^, there exists an integer L{k) in O(logfc) such that either no 
column is fired on the right of L{k), or columns L(k) and L(k) + 1 both are fired 
(and therefore, the property of the second item above applies). 

Informally, that means that we have the emergence of a regular behavior, after a short 
transitional and complex phase. 

These results give a better understanding of avalanches for sufficiently large columns. 
We hope that in future work, they will help us in the approach of the structures of fixed 
points TT^k). 

1.3 The context 

The problem of describing and proving regularity properties, experimentally checked, for 
models issued from basic dynamics is really a present challenge for physicists, mathemati- 
cians, and computer scientists. There exists a lot of conjectures, issued from simulations, 
on discrete dynamical systems with simple local rules (sandpile model [3] or chip firing 
games, but also rotor router [13], the famous Langton ant [7] [8]...) but very few results have 
actually been proved. As regards KSPM(Z)), the prediction problem (namely, the problem 
of computing the fixed point 7r(A;) knowing 7r(fc — 1)) has been proven in [l4j to be in NC^ 
for the one dimensional case (the model of our purpose), which means that the time needed 
to compute an avalanche is in 0(log'^ N) where A'^ is the number of grains, and P-complete 
when the dimension is > 3. A recent study ([lOj) showed that in the two dimensional 
case the avalanche problem (given a configuration a and a column i on which we add one 
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grain, does it have an influence an index j?) is P-complete, which points out a inherently 
sequential behavior. 

This study will provide tools to understand sand pile evolution. We hope that those 
tools form a basis to obtain some good descriptions of fixed points ir{k), but are also 
deeply related with other subjects around sand piles such as unit elements of abelian group 
structures presented in [2] and [4]. 

2 Avalanche process in the general case 

This section begins with a first glance at avalanches, allowing notation simplifications. 
Then avalanches are studied in details, leading to a simplified description of its behavior. 

Proposition 1. For each strategy s such that tt{N)^'^ A 7r(A^ + 1) and each i G N , we 
have \s\i £ {0, 1}. 

Proof. Let s = (si, . . . , st) be a strategy such that 7r(A^)-'''' A 7r(A^ + 1). We have to prove 
that, for 1 < / < m < r, we have si ^ Sm (obviously, \s\i > for all i). To do it, we prove 
by induction t <T that for 1 < / < m < t, we have si ^ Sm- 

For initialization this is obviously true for t = 1. Now assume that the condition is satisfied 
for an integer t such t < T, and let z be a column such that there exists an integer / < t 
such that i = si- Let a be the configuration such that 'k{N)^^ -4 . . . A cj. 
Notice that the transitions which can possibly change the value of the current configuration 
at i could be: i (which decreases the value by D units), i + 1 (which increases the value by 
D — \ units) or i — D + 1 (which increases the value by 1 unit). 

Thus we have ai < 7r{N)f - D + D -l + l since by definition, between 7r(A'^) and a, exactly 
one transition has occurred in i, at most one transition has occurred in z + 1, and at most 
one transition has occurred in i — + 1. For i > 1, we get < Tr{N)i. On the other hand, 
since 7r(A^) is a fixed point, we have: vr(A^)j < D, which guarantees that st+i ^ i- For i = 0, 
there is no possible transition in i — D + 1, thus we get do < TT{N)jl^ — D + D — 1, which 
is do < 7r{N)Q + 1 — D + D — 1. Thus ctq < ■7r(A^)o < D which also gives: st+i / 0. 
This ensures that the result is true for t + 1, and, by induction, for T. 

When talking about an avalanche s, lemma [1] allows us to write z E s instead of \s\i = 1 
without lose of information. We denote by S[„ „] the subsequence of s from u to f included. 

We will now study avalanches in details. For D = 2, i.e. the classical SPM, avalanches 
are quite simple, the added grain moves rightwards until it finds a stable position. For 
D > 2, the situation is more complex, and needs a precise study, given by the following 
lemma. 

Lemma 1. Let s = {si, . . . , st^.) be the fc*^ avalanche. Let rt = maa;{sf',t' < t}. 

— Assume that st+i < rt. Then St+i is the largest column number i such that i < rt and 
i ^ S[i,t]- Moreover, we have: rt — St+i < D — 1. 
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— Assume that st+i > Vf. We have: st+i — rt < D — 1. 

Proof. We first order fired columns by causality. Precisely, a column i has two potential 
predecessors, which are i + 1 and i — D + 1. State i = Su- These columns are really 
predecessors of i if they are elements of i-6 if they are fired before i. By this way, 

using the transitive closure, we define a partial order relation (denoted by <causX ™ fired 
columns for s. 

Now, consider the set At+i of ancestors of st+i (i.e. the set of columns i such that 
i <caus. st+i) and, the set St of which have st as common ancestor (i.e. columns i such 
that St <caus. i)- We necessarily have rt G ^t+i. Otherwise, we have ^t+i Pi St = ^, and 
this allows another strategy s' , constructed from s postponing the transitions at rt and 
elements of St after the transition on st+i- This contradicts the fact that s is leftmost. 

Let {io,ii, ...ip) be a finite sequence such that iq = rt, ip = st+i and, for each j with 
< j < p, ij is a predecessor of ij+i- Such a sequence exists since rt G ^t+i- One easily 
proves by induction that ij = rt — j: this is true for i = 0. Assume it is true until the 
integer j < p. We have either ij^i = ij — 1 or ij^i = ij + D — 1. But from the induction 
hypothesis, ij + D — 1 is an ancestor of ij^i or has not yet been fired, thus ij+i = ij — 1. 
This gives that s^+i is the largest column number i such that i < rt and i ^ s^i^t]- 

Now if we assume, by contradiction, that p > D — 1, then rt — D + 1 is not a predecessor 
of rt, which yields that rt has no predecessor, which is a contradiction. This gives the 
inequality of the first item. The second item is obvious, since st+i ha a unique predecessor 
which is st+i — D + 1. 

Lemma [U induces a partition of fired columns between those which make a progress {i.e. 
increases the greatest fired column) and those which do not. This distinction is important 
in further development, so let us give progress firings a name. Let s = (si, . . . , st) be an 
avalanche, a column st is called a peak if and only if st > maxs^i^t-i]- 

Remark 1. Two peaks p ^ q can be compared using chronological (<t) or spatial {<s) 
orders. Nevertheless, by definition of peaks we obviously have p <t q <=^ p <s q- 

Lemma 2. Let s be the k^^ avalanche. Assume that there exists a column I, such that for 
each column i with l<i<l + D — l,i£s, and a column i' such that i' > I + D — 1 and 
i' £ s. Let I' be the lowest peak such that V >l + D — 1. 

There exists a time t such that: 

( for all i with I' — D + l<i<l', i £ S[i 

\for all i with I' < i, i ^ s^i^t] 

Moreover I' is the lowest integer such that I' > I + D — 1 and aj, = D — 1. 

Proof. Let to be the time when st^ = I', i.e. the first time such that st^ > I + D — 1, and 
let j be the largest integer such that, for < j' < j, we have stg+j' = st,, — j' ■ Let us state 
* = *o + j- We have j < D - I. 
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Let fj* denote the configuration obtained from ir{k — 1) via Let i, witli i < I', sucli 
tfiat i ^ S[i.t]- We claim that we have : i ^ s. To prove it, we prove by induction that for 
any t' > t, i ^ ^[i,t']- Assume that this is satisfied for a fixed t'. This means that all the 
transitions of S[(_|_i^j'] are done on columns larger than Thus, aj = aj and no transition 
is possible on i for a* since s is leftmost (the only potential column to be fired is Sfp — j — 1, 
but by assumption, either this column has been previously fired, or it cannot be fired by 
definition of j, according to lemma [1]) . 

By contraposition, it follows that for each column i with / <i</ + D — l,we have i G ^j. 
A simple (reverse sense) induction shows that, for l + D — l<i<l' we have i G s, since 
by hypothesis i + 1, and i + 1 — D both are in s. Thus, by contraposition of the claim 
above, for l + D — l<i<l', we have i S -s^i^i]. This gives the the fact that for all i with 
l'-D + l<i<l',i£ s[i^t] 

The fact that for all i with I' < i, i ^ ^['^,t] is trivial, by definition of to and t. 
We have /' > 1 + D — 2 and aj, = D — 1. assume that there exists /" < /' satisfying the same 
properties. Notice that for to < t' < t we have s'f- > I' — D — 1. Thus the time ti such that 
st-^ = I" — D + 1 is such that ti < to- That means that /" should have been fired before to; 
a contradiction. 

Lemma [2] describes in a very simple way the behavior of avalanches. Thank to it, the 
study of an avalanche can be turned into a pseudo linear execution, in which transitions 
are organized in a clear fashion: 

Theorem 1. Let s = (si, . . . , st) be the k^^ avalanche and {pi ■ ■ ■ ,Pq) be its sequence of 
peaks. Assume that there exists a column I, such that for each column i with I < i < l+D—1, 
i £ s. Then for any column p such that p > I + D — 1, 

p is a peak of s <^=^ p ^ Pq + D — I and 7r{k — l)p = D — 1 

Furthermore, Let Pi = St, with pi >l + D — \, be a peak. Then 

T >t + Pi — Pi-i — 1 and for all t' s.t. t < t' < t + pi — Pi-i — 1, = St'-i — 1 

A graphical representation of this statement is given on figure [3l 

Proof. The first part is a straight induction on lemma [2j 

The second part follows an induction summed up in the following fact: any column i such 
that 7r(/c — l)i < D — 1 must wait for its right neighbor i + 1 to be fired, and it should be 
fired when both i + 1 and i — D + 1 has been fired (besides, i — D + 1 has already been 
fired). Since any of such i is fired to reach a fixed point, T >t + pi — pi-i — 1. 

In the light of remark [H theorem [1] easily allows us to compute the right part of k^^ 
avalanche (from column I + D — 1), only knowing 7r(A; — 1). The sequence of peaks is 
computed as follows. The first one is the lowest column i greater or equal to I + D — \ 
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Fig. 3: Illustration of Theorem [T] with D = 6. Surrounded columns I to I + D — 2 are 
supposed to be fired. Black column is the greatest peak strictly lower than L + D — 1. A 
column is grey if and only if its value is -D — 1. Following arrows depicts the avalanche. 



such that Tr{k — l)i = D — 1. Then, given a peak i, the next one is the lowest j such that 
iT{k — l)j = D — 1 and j — i < D — 1. If such a j does not exist, then there is no more peak 
and i is the largest fired column. 

We can distinguish two movements within an avalanche: before a certain column it has 
an unknown behavior, and from that column to the end the behavior is pseudo local, in 
the sense that when an index is fired ahead (on the right) then any 'hole' is filled before 
the progress can continue. 

An important direct implication of theorem [T] is that if there exists a column I such 
that for the k^^ avalanche s^, we have for all l<i<l-\-D — l,i£s, then for all j such that 
l + D — l<j< maxs^, we have j — {D — + 1 G s'^ and therefore 7T{k)j = 7r{k — 

Intuitively, this equality hints some similarity between successive avalanches. 

Note that previous results also apply for a grain addition on column of any fixed 
point configuration of KSPM(L'). 

This study constitutes a simplified understanding of the behavior of avalanches, which 
we hope will be helpful toward the description of fixed points. As motivated above, next 
subsection studies, for KSPM(3), the previous result hypothesis that for an avalanche s 
there exists a column I such that for all I < i < I + D — 1, i is element of s. 



3 Short transitional phase when D = 3 

In this section we prove that in KSPM(3), there exists a column 1{N) in 0(log N) such that 
lemmas hypothesis is verified for any avalanche s^, with k < N, such that maxs'^ > l{k). In 
other words, considering the first avalanches, from a logarithmic column, we can apply 
theorem [T] and consider avalanches pseudo locally, as described on figure [31 Here is the 
statement: 

Proposition 2. Let s be the k^^ avalanche of KSPM(3). There exists a column l{k) in 
0{log k) such that for any k, when max{j|j G s} > l{k), l{k) and l{k) + 1 both are 
elements of s. 

Proof. Let i be a fixed column. If i,z + 1 S , then, theorem [T] states that for all i' such 
that i < i' < max{j|j E s}, we have i',i' + 1 G s. If j, j + 1 ^ s, then, from Proposition [H 
we have max{j|j G s^} < i. 
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Let j be a fixed positive integer. Assume that the avalanche s fires 2j but not 2j — 1. 
Prom remarks above, columns 0, 2, 4, 2j are fired in s while columns 1, 3, 5, 2j — 1 are 
not fired. 

If a column i + 1 is fired while i is not, then we necessarily have Tr{k — l)i = 0, since the 
firing in z + 1 increases the value in column i from 2 units. Moreover, if the column i + 1 
is fired while i + 2 is not, then we necessarily have Tr{k — l)i+i = 2, since the i + 1 receives 
at most one grain, by preceding firings. 

On the other hand, obviously, the assumption on j enforces that Tr{k — l)o = 2 This 
yields that (2, 0)-'^^ is a prefix of iT{k — 1). 

We have the following fact : 

Fact: There exists constant numbers A and B, with ^ > such that if a configuration 
7r(iV) has a prefix of the form (2, 0)^' then N > AA^ + B 

This is obtained by the linear algebra analysis below. This gives the result of the 
proposition. 

Let 7r(iV) = ((Jo, cxi, . . . ) be the configuration and a = (oq, ai, . . . ) be its shot vector i.e. 
the sequence a = (ao, ai, 02, . . . ) where is the number of times the column i has be fired 
in the N first avalanches. According to the iteration rule we have the relation: 



= 0-1-2 — 3ai + 2aj+i 

i.e. 

a.i+1 = o.i-2 + 3aj) 

/ 1 \ 

We state ^ = 1. We denote by Vi the column vector such that and vj = 

\-l/2 3/2/ 

(0, 0, (Ti/2), and Ui the column vector such that uf = (ai_2, aj-i, a^) (with the convention 
that is the row vector obtained by transposition of the column vector u) . The equality 
above can be algebraically written in 

Ui+i = Aui + Vi 

By iteration we get : 

iii+2 = A?Ui + Avi + Vi+i 

What we want is a2i = 2 and (T2i+i = so = (0, 0, 1) and 'u^+i = (0, 0, 0). With this 
specification, we get : 

«2(i+l) = ^^«2i + b 
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with = (0, 1,3/2). From this last relation we will deduce a condition on N to get 
the sequence (2,0)-' . 



Let us first find a new basis to get the matrix A on Jordan canonical form. The char- 
acteristic polynomial of A is ^{2x + l)(x — 1)^ and its eigenvalues are —5 of algebraic 
multiplicity 1 and 1 of algebraic multiplicity 2. Since dim(ker(^ — Id)) = 1 the Jordan 
canonical form of A is 

/II 

^Jordan =01 

V0 0-l/2y 

And the new basis E' = (e'^, 63, 63) according to the canonical one E = {ej = (1, 0, 0), 62 = 
(0, 1, 0), 63 = (0, 0, 1)) is given by the linear relations: 

e[ = 61 + 62 + 63 

62 = 62 + 263 

63 = 461 - 262 + 63 

Let p be the linear mapping consisting in the projection on the line D3 generated by 63 
according to the direction the plane Pi, 2 generated by e'^ and 63. For any vector u, we have: 
p{A'^u) = p{A'^{p{u) + u—p{u)) = p{A'^ {p{u)) + p{A'^ {u — p{u)) . Notice that, by definition of 
p, u — p{u) is element of Pi,2, and, therefore A{u—p{u)) and A^{u—p{u)) also are elements 
of Pi, 2- This yields that p{A^{u — p{u))) is null. On the other hand, since p{u) is element of 
1)3, A{p{u)) = ^p{u) and A^{p{u)) = \p{u)\ thus p{A^{p{u)) = \p{u). As a conclusion, 
we get p{A^u) = \p{u), which, in particular, allows the following equalities: 

^(""2(1+1)) =p{A^U2i + b) 

= p{A^U2i) + p{h) 
= \p{u2i) - i^e^j 

Let V be the unique vector collinear with 63 satisfying the equation 

1 1 / 

V = -V 60 

4 18 ^ 

\. e. V = ^63. Remember that p{v) = v. We have 

P(«2(i+1) -V)= p{u2(i+i)) - V 

= lpiu2i) - js^s - {\v - ^e'^)) 

= \p{u2i - V) 

This gives by induction: 

p{uo - v) = 4^p{u2j - v) 
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Now we specify the sequence of vectors Ui, assuming that values Oj are the shot vectors 
of a configuration a beginning by (2,0)-^. (For convention we also state a_2 = N and 
a_i = 0, thus we have Uq = {N,0,aQ), w[ = (0,ao,ai) and uJ = (aj-_2, Oi-i, for i > 2). 
An easy computation gives that: p{uo — v) = ^+°°+27 
Let Xj be defined by p{u2j — v) = xje'^- We obtain the equality : 



27 - 



J 



Obviously ao < ^ = -j, which ensures that N + ao + < ^ + 1. 



Furthermore, we necessarily have xj > 0, and, from lemma [3] proved on the bounce, 
each element of p(Z^) is a multiple of ce'^, where c is a positive constant. If v is element 
of p{Z^) we can conclude that Xj > c. If v is not element of p{T?), we can conclude that 
Xj > min{\ck + ^|, /c G Z}. In any case, there exists a positive real d, not depending on j, 
such that Xj > d. 

We conclude that ^ + 1 > A^d, which gives > 4-' — | to get a sand pile of the 
form a = (2,0)^0-'. 

We now give the lemma used in the previous proof. 

Lemma 3. (constant steps) There exists apositive realc such that p^^i^iW?) = {ice'^^i G Z}. 

Proof. The set of reals r such that there exists an element x in such that Pe>,^{x) = re!-^ 
is obviously a group. So we only have to prove that this group is discrete, i.e. that there is 
no sequence (r„)„gz of positive reals such that lim (r„) = 0. 

Assume, by contradiction, the existence of such a sequence, and let (x,i)„6z be a sequence 
of vectors such that, for each integer n, Pe'JyXn) = rn. 

A key-point is that vectors e'^, ^"^^ ^3 have integer components, so we can state x„ = 
ane'i+^nCa+Cneg. The sequence (x'„)ngz defined by x'^ = (a„- [anJ)e'i + (6„- [6„J)e2+c„e3 
also is a sequence of integer vectors such that for each integer n, Pe'^i^'n) = rn. Moreover 
this sequence is bounded. Thus (x^)^ S Z takes a finite number of values, which enforces 
that the sequence (r„)„ G Z also takes a finite number of values, which is a contradiction. 

For KSPM(3), after a short transitional of logarithmic length, hypotheses of theorem 
[1] are verified , and the study of avalanches can be turned into a pseudo linear process. 
Note that a trivial framing of the maximal non-empty column e(A^) of a fixed point with 
A^ grains shows that e(A^) is in fi{y/N). As a consequence, pseudo local process stands for 
the asymptotically complete behavior of avalanches. 

Unfortunately, the approach above does not hold for D > 3. The main reason is that, 
for D = 2> unfired columns induce a very particular and periodic prefix ((2, 0)-') on config- 
urations. From Z) = 4, the structure of such a possible prefix is more complex and we did 
not yet get a tractable characterization of those prefixes. 
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4 Perspectives 



In this paper we described avalanches as pseudo local processes from a certain column I. 

We proved this column to be logarithmic in the number of grains N for KSPM(3), 
leading to an asymptotically complete description of avalanches in that case. Simulations 
for other parameter D suggests that the same outcome also holds. 

The pseudo local process description involves some properties on avalanches, which we 
hope will be useful toward the study of fixed points shape. For an avalanche s, a particularly 
interesting consequence is that two STicccssivc fixed points are equal from I + D — 1 to 
(max .s) — 1, which hints that next avalanche reaching this part of the configuration may have 
a similar behavior. This would lead to a knowledge on the likeness of successive avalanches 
and therefore a foresee on the shapes of fixed points. Further work may concentrate on 
this point, where the main purpose is to go ahead iterating evolution rules, and to describe 
fixed points with a plain formula. 

References 

1. Per Bak, Chao Tang, and Kurt Wioscnfcld. Solf-organizod criticality. Phys. Rev. A, 38(1) :364 374, Jul 
1988. 

2. Michael Creutz. Cellular automata and self organized criticality. In in Some New Directions in Science 
on Computers, 1996. 

3. Arnaud Dartois and Clemence Magnien. Results and conjectures on the sandpile identity on a lattice. In 
Michel Morvan and Eric Remila, editors, Discrete Models for Complex Systems, DMCS'03, volume AB 
of DMTCS Proceedings, pages 89-102. Discrete Mathematics and Theoretical Computer Science, 2003. 

4. Deepak Dhar. Self-organized critical state of sandpile automaton models. Phys. Rev. Lett., 64(14);1613- 
1616, Apr 1990. 

5. Jerome Olivier Durand-Lose. Parallel transient time of one-dimensional sand pile. Theor. Comput. 
Sci., 205(1-2):183-193, 1998. 

6. Enrico Formenti, Benoit Masson, and Theophilos Pisokas. Advances in symmetric sandpiles. Fundam. 
Inform., 76(1-2):91-112, 2007. 

7. Anahi Gajardo, Andres Moreira, and Eric Goles. Complexity of langton's ant. Discrete Applied Math- 
ematics, 117(l-3):41-50, 2002. 

8. David Gale, Jim Propp, Scott Sutherland, and Serge Troubetzkoy. Further travels with my ant. Math- 
ematical Entertainments column, Mathematical Intelligencer, 17:48-56, 1995. 

9. Eric Goles and Marcos A. Kiwi. Games on line graphs and sand piles. Theor. Comput. Sci., 115(2):321- 
349, 1993. 

10. Eric Goles and Bruno Martin. Computational Complexity of Avalanches in the Kadanoff Two- 
dimensional Sandpile Model. In TUGS, editor. Proceedings of J AC 2010 Joumees Automates Cellulaires 
2010, pages 121-132, Turku Finland, 12 2010. F.1.1. 

11. Eric Goles, Michel Morvan, and Ha Duong Phan. The structure of a linear chip firing game and related 
models. Theor. Comput. Sci., 270(l-2):827-841, 2002. 

12. Leo P. Kadanoff, Sidney R. Nagel, Lei Wu, and Su-rnin Zhou. Scaling and universality in avalanches. 
Phys. Rev. A, 39(12) :6524-6537, Jun 1989. 

13. Lionel Levine and Yuval Peres. Spherical asymptotics for the rotor-router model in z d. Indiana Univ. 
Math. J, pages 431-450, 2008. 

14. Cristophor Moore and Martin Nilsson. The computational complexity of sandpiles. Journal of Statistical 
Physics, 96:205-224, 1999. 10.1023/A:1004524500416. 

15. Thi Ha Duong Phan. Two sided sand piles model and unimodal sequences. ITA, 42(3):631-646, 2008. 



12 



